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Abstract. We study various aspects of the dynamics induced by integer matrices on the 
invariant rational lattices of the torus in dimension 2 and greater. Firstly, we investigate 
the orbit structure when the toral endomorphism is not invertible on the lattice, character- 
ising the pretails of eventually periodic orbits. Next we study the nature of the symmetries 
and reversing symmetries of toral automorphisms on a given lattice, which has particular 
relevance to (quantum) cat maps. 



1. Introduction 

Toral automorphisms or cat maps, by which we mean the action of matrices M € GL(d, Z) 
on the d-torus T rf , are a widely used and versatile class of dynamical systems, see |45[ [28] for 
some classic results in the context of ergodic theory. Of particular interest are the hyperbolic 
and quasi-hyperbolic ones, which are characterised by having no root of unity among their 
eigenvalues. All periodic orbits of such automorphisms lie on the rational (or finite) invariant 
lattices L n = {x £ T d | nx = mod 1}, which are also known as the n-division points. One 
can encode the possible periods of a toral automorphism M on T d via the dynamical zeta 
function in a systematic way, which is always a rational function [HJ [18] . The literature on 
classifying periodic orbits of toral automorphisms when d = 2 is vast (compare |21 [ I25 [ [37] and 
references therein). An extension beyond d = 2 is difficult due to the fact that the conjugacy 
problem between integer matrices is then much harder because (unlike d = 2) no complete 
set of conjugacy invariants mod n is known. Therefore, our focus will also be on d = 2, with 
occasional extensions to higher dimensions. 

The larger ring Mat (ci, Z) of toral endomorphisms (which includes integer matrices without 
integer inverses) has received far less attention [2] [14"1 [8] , particularly those in the complement 
of GL(d, Z). Note that the resulting dynamics induced by M G Mat(d, Z) \ GL(d, Z) on a 
finite lattice L n may or may not be invertible. In the latter case, beyond periodic orbits, there 
exist eventually periodic orbits which possess points that lead into a periodic orbit. We call 
these points and the periodic point to which they attach the 'pretails' to the periodic orbit 
(see Eq. (I15p for a formal definition). The action of M induces a directed graph on L n (e.g. 
see our three figures below). Alternatively, the pretails can be combined to form a rooted tree 
which is a characteristic attribute to any pair (M, L n ). 

As well as their interest from a mathematical viewpoint, toral automorphisms also have 

been well-studied from a physics perspective, in particular as quantum cat maps (see [30, 
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\19\ [32] and references therein). Here, the action of the integer matrix on a rational lattice 
L n , for some n, is all-important as quantum cat maps and their perturbations are built from 
(classical) cat maps and their perturbations restricted to a particular rational lattice (called 
the Wigner lattice in this instance). There has been recent interest in dealing with so-called 
pseudo-symmetries of quantum cat maps that are manifestations of local symmetries of cat 
maps restricted to some rational lattice [301 US E2] • Although, in the context of quantisation, 
matrices from the group Sp(2d, Z) play the key role, we prefer to work with the larger group 
of unimodular integer matrices and consider the former as a special case. 

The main aims of this paper are twofold: (i) to elucidate the orbit structure of toral en- 
domorphisms on rational lattices, equivalently the periodic orbits together with the related 
pretail tree structure; (ii) to further characterise the nature of symmetries or (time) reversing 
symmetries of toral automorphisms, these being automorphisms of the torus (or of a rational 
lattice) that commute with the cat map, respectively conjugate it into its inverse. 

We expand a little on our results, where we refer to the actual formulation below in the 
paper. The results are readable without the surrounding notational details. 

With respect to aim (i), Section [3] characterises the splitting of L n into periodic and even- 
tually periodic points under a toral endomorphism M. Every periodic point has a pretail 
graph isomorphic to that of the fixed point (Corollary [T]), which is trivial if and only if M is 
invertible on L n . In general, the pretail tree codes important information on the action of M. 
One question in this context is whether all maximal pretails have the same length, for which 
we give a partial answer via a sufficient condition on ker(M) in Proposition [3l Given M, the 
lattice L n can be decomposed into into 2 invariant submodules, one of which captures the 
invertible part of M and the other the nilpotent part. This way, we are able to characterise 
the dynamics that is induced by M on L n in the case of n = p r , p prime, in Corollary [3] and 
Lemma [3j 

Our contribution towards aim (ii) continues the investigations from [9J [TOl [T2] . The key 
quantity for integer matrices of dimension 2 is the mgcd (see Eq. 1131 below), and one conse- 
quence of [121 Thm. 2] is that M G SL(2, Z) is always conjugate to its inverse on L n , for each 
n G N. The conjugating element - called a reversing symmetry or reversor - is an integer 
matrix that has an integer matrix inverse on L n , which typically depends on n. In this way, 
any SL(2,Z) matrix that fails to be conjugate to its inverse on the torus (e.g. M = (| 1 g) 
from [£l Ex. 2]) is still conjugate to its inverse on every rational lattice. In [T2|, we did not 
consider the nature of the reversor on the lattice. Theorem [T] of Section [4] establishes that 
it is an orientation-reversing involution, what is called an anticanonical (time-reversal) sym- 
metry in the language of [30]. Section [4.21 uses normal forms of GL(2,F p ) to characterise the 
symmetries and possible reversing symmetries of such matrices; the underlying structure of 
the conjugacy classes of GL(2,F p ) is summarised in Table l4~2l The symmetry structure has 
some extensions to higher dimensions (Section 14.3ft and to general modulus n (Section 14.4ft . 
Section 14.51 presents some results for the case when M G GL(d, F p ) has a root in the same 
group. 
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The structure of the paper is as follows. Section [2] summarises some properties of integer 
matrices that we use later in the paper, with some reformulations or slight generalisations 
that we find useful. In particular, throughout the paper, we formulate the results for arbitrary 
dimension whenever it is possible without extra complications, though this is not our main 
focus. As described, Section addresses aim (i) above, while Section H] deals with aim (ii). In 
the Appendix, we briefly discuss two classic examples of toral automorphisms for d = 2 and 
some aspects of their dynamics. 

2. Preliminaries and powers of integer matrices 

The purpose of this section is to summarise important properties of and around integer 
matrices that are needed later on, with focus on those that are not standard textbook material. 
At the same time, we introduce our notation. For general background on integer matrices 
and their connections to algebraic number theory, we refer to the classic text by Taussky |43| . 

2.1. Lattices, rings and groups. The most important lattices on the torus T d = M d /Z d , 
which is a compact Abelian group, consist of the n-division points 

(1) L n := {x G T d j nx = (mod 1)} = {(§■,. . . , | < h < n for all 1 < i < d] , 

with n G N. Clearly, the L n are invariant under toral endomorphisms (with the action of 
the representing matrices taken mod 1). It is sometimes easier to replace L n by the set 
L n := {(fci, . . . , kd) | < ki < n}, with the equivalent action of M defined mod n. This also 
applies to various theoretical arguments involving modular arithmetic. Consequently, we use 
L n (with action of M mod 1) and L n (with action mod n) in parallel. 

Our discussion will revolve around the residue class ring Z/nZ with n G N, which is a 
principal ideal ring, but not a domain, unless n = p is a prime. In the latter case, Z/pZ = ¥ p 
is the finite field with p elements, while the ring has zero divisors otherwise. For general n, 
the unit group 

(Z/nZ) x = {1 < m < n \ gcd(m,n) = 1} 

is an Abelian group (under multiplication) of order cp(n), where <j> is Euler's totient function 
from elementary number theory |26j. In general, it is not a cyclic group. 

The integer matrices mod n form the finite ring Mat(d, Z/nZ) of order n d2 . The invertible 
elements in it form the group GL(d,Z/raZ) = {M G Mat(d,Z/nZ) | det(M) G (Z/nZ) x }. If 
n = p^ 1 ■ ■ -pf is the standard prime decomposition, one finds 

(2) |GL(d,Z/nZ) 
where 

(3) |GL(d,F p )| = (p d 

is well-known from the standard literature [331134], Formula (J2J follows from the corresponding 
one for n = p r via the Chinese remainder theorem, while the simpler prime power case is a 
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consequence of the observation that each element of a non-singular matrix M over Z/p s Z can 
be covered (independently of all other matrix elements) by p elements in Z/p s+1 Z without 
affecting its non-singularity. 

Let us finally mention that SL(n, Z/nZ), the subgroup of matrices with determinant 1, is 
a normal subgroup (it is the kernel of det: GL(n,Z/raZ) — > (Z/nZ) x ). The factor group is 

GL(n, Z/nZ) / SL(re, Z/nZ) ~ (Z/nZ) x 

and thus has order 4>{n). 

2.2. Orbit counts and generating functions. The orbit statistics of the action of a matrix 
M 6 Mat(d, Z) on the lattice L n is encapsulated in the polynomial 

(4) Z n (t) = \{{l-t m )^\ 

meN 

(n) 

where c m denotes the number of m-cycles of M on L n . Recall that, if a m and c m denote the 
fixed point and orbit count numbers of M (dropping the upper index for a moment), they 
are related by 

(5) a m = ^dc d and c m = — X^(lf) a ^ 

777 

d\m d\m 

where n(k) is the Mobius function from elementary number theory [26j . 

Despite the way it is written, Z n is a finite product and defines a polynomial of degree at 
most n d . Note that the degree of Z n can be smaller than n d (as the matrix M need not be 
invertible on L n ), but Z n (t) is always divisible by (1— t), because is a fixed point of every M. 
The polynomials Z n are closely related |12[ [36] to the zeta function of toral endomorphisms, 
which can be calculated systematically; compare [8] and references therein. Dynamical zeta 
functions give access to the distribution and various asymptotic properties of periodic orbits 
[El [40], an d also relate to topological questions; compare [23] for a systematic exposition of 
the latter aspect in a more general setting. Further aspects on the asymptotic distribution of 
orbit lengths on prime lattices can be found in [29] . 

2.3. Matrix order on lattices and plateau phenomenon. Assume that M is invertible 
on L n (hence also on L n ). Then, its order is given by 

(6) ord(M,n) := gcd{m G N | M m = 1 modn}. 

Clearly, ord(M, 1) = 1 in this setting. When M is not invertible on L n , the definition results 
in ord(M, n) = 0; otherwise, ord(M, n) is the smallest m € N with M m = 1 mod n. 

Let M G GL(d, Z) be arbitrary, but fixed. To determine ord(M, n) for all n > 2, it suffices 
to do so for n an arbitrary prime power, since the Chinese remainder theorem [26j gives 

(7) ord(M,n) = lcm(ord(M,Pi 1 ), . . . , ord(M,j$)) 

when n = p r ^ ■ ■ -pf is the prime decomposition of n. It is clear that oid(M,p r )\ ord(M, p r+l ) 
for all r G N, see also [141 Lemma 5.2]. 
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Let us now assume that M G Mat(d, Z) is not of finite order, meaning that M k ^ 1 for all 
k G N, which excludes the finite order elements of GL(d, Z). If p is a prime, we then obtain 
the unique representation 

(8) M ord ( M ' p ) = 1 + p s B 

with s G N and an integer matrix B ^ mod p. Starting from this representation, an 
application of the binomial theorem for powers of it, in conjunction with the properties of the 
binomial coefficients mod p, gives the following well-known result. 

Proposition 1. Let M £ Mat(d, Z) be a matrix that is not of finite order. Fix a prime p 
that does not divide det(M), and let s be defined as in Eq. ([8]). 

When p is odd or when s > 2, one has ord(M,p l ) = oid(M,p) for 1 < i < s, together with 
ord(M,p s+i ) = p i ord(M,p s ) for all i G N. 

In the remaining case, p = 2 and s = 1, one either has ord(M, 2 r ) = 2 r ~~ 1 ord(M, 2) for all 
r G N, or there is an integer t > 2 so that ord(M, 2 l ) = 2ord(M, 2) for 2 < i < t together 
with ord(M, 2 t+i ) = 2* ord(M, 4) for all i G N. □ 

In what follows, we will refer to the structure described in Proposition [1] as the plateau 
phenomenon. Such a plateau can be absent (p odd with s = 1, or the first case for p = 2), it 
can be at the beginning {p odd with s > 2), or it can occur after one step (p = 2 when t > 2 
exists as described), but it cannot occur later on. 

Proposition [1] is a reformulation of |144 Thms. 5.3 and 5.4], which are originally stated for 
M G GL(2, Z). As one can easily check, the proofs do not depend on the dimension. Similar 
versions or special cases were also given in [13] and [H] (with focus on SL(2, Z)-matrices), in 
|37j (for the order of algebraic integers), in [44J (for the Fibonacci sequence), in [15] (for linear 
quadratic recursions) and in [22] and [16] (for general linear recursions). Let us also mention 
that, based on the generalised Riemann hypothesis, Kurlberg has determined a lower bound 
on the order of unimodular matrices mod N for a density 1 subset of integers iV in |31j . 

2.4. Powers of integer matrices. Consider a matrix M G Mat(<i, Z) with d > 2 and 
characteristic polynomial P M (z) = det(zl — M), which (following |46| ) we write as 

Pm( z ) = z<1 — c x z d 1 — c 2 z d 2 — ... — c d _ 1 z — c d , 

so that c d = (— l) d+1 det(M). Let us define a recursion by u = = . . . = u d _ 2 = and 
Uj x = 1 together with 

d 

(9) u m = ^2 c i u m _ i = c 1 tt m _ 1 + c 2 u m _ 2 + . . . + c d u m _ d 

1=1 

for m > d. This results in an integer sequence. Moreover, when c d ^ 0, we also define 

U m = c d 1 ( U m+d ~ C l U m+d-l ~ ■ ■ ■ ~ c d-l u m+l) 
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for m < — 1. In particular, since d > 2, one always has u_ 1 = l/c d and u_ 2 = — c rf _ 1( /c^, 
while the explicit form of u m with m < —2 depends on d. Note that the coefficients with 
negative index are rational numbers in general, unless \c d \ = 1. 

The Cayley-Hamilton theorem together with ([9]) can be used to write down an explicit 
expansion of powers of the matrix M in terms of M k with < k < d — 1, 

d-i 

(10) M m = Y,lt )M ^ 

where the coefficients satisfy 7^ = 5 m e (for < £, m < d — 1) together with the recursion 

(11) i { r 1] =c d -n d n \ + itX 

for n > d — 1 and < ^ < d — 1, where 71^ := 0. In particular, = c d _ e . The coefficients 
are explicitly given as 

I d-e-i 

(12) li = '^2 c d-i u m-£-l+i = u m+d-l-l ~ ^ C d-l-i u m-l+i ) 

i=0 i=l 

where m > d and the second expression follows from the first by Q . Formulas (|10p and (|12p 
can be proved by induction from M d = c x M d ~ x + c 2 M d ~ 2 + . . . c d _ x M + c d l. Eq. dTUJ) holds 
for all m > in this formulation. 

When det(M) 7^ 0, the representation (|12p also holds for m < (i, as follows from checking 
the cases < m < d together with a separate induction argument for m < 0. In particular, 
one then has 

M _1 = c d n_ 2 l + (c d _ 1 «_ 2 + c d n_ 3 )M + (c d _ 2 n_ 2 + c d _!«_3 + c d u_ 4 )M 2 
+ . . . + (c 2 u_ 2 + c 3 u_ 3 + . . . + c d u_ d )M d - 2 + u„ 1 M d - 1 , 

which is again an integer matrix when \c d \ = 1. 

2.5. Results for d = 2. Let us look at matrices from Mat(2,Z) more closely, and derive one 
important result by elementary means. Consider M = ( a c b d ), set D := det(M), T := tr(M) 
and define the matrix gcd (or mgcd for short) as 

(13) mgcd(M) := gcd(6, c,d — a), 

which is another invariant under GL(2,Z) conjugation. Its special role becomes clear from 
the following result, which is a reformulation of |12} Lemma 2 and Thm. 2]. This will lead to 
Corollary [2] below. 

Lemma 1. Two matrices M,M' G Mat(2,Z) that are GL(2, r L)-conjugate possess the same 
mgcd, as defined in Eq. (|13|) . More generally, the reductions modulo n of M and M' are 
GL{2, r L/nL)-conjugate for all n > 2 if and only if the two matrices share the same trace, 
determinant and mgcd. □ 



V 

Returning to matrix powers, formula (jlOp simplifies to 
(14) M m = u m M-Du m _ x t, 

where now it = 0, u x = 1 and u m+1 = Tu m — Du m _ 1 for m G N; see [I2j Sec. 2.3] for details. 
Let n 6 N and assume gcd(n, D) = 1. This allows us to introduce 

n(n) := period of (u m ) m >o mod n 

which is well-defined, as the sequence mod n is then indeed periodic without 'pretaiF. Recall 
that (it m )m>o mod n must be periodic from a certain index on, as a result of Dirichlet's 
pigeon hole principle. Since D is a unit in Z/nZ, the recursion (|14p can be reversed, and 
(u m )m>o mod n must thus be periodic, with k(h) being the smallest positive integer k such 
that u k = and u k+1 = 1 mod n. 

One can now relate n(n) and ord(M, n) as follows, which provides an efficient way to 
calculate ord(M, n). 

Proposition 2. Let M E Mat(2,Z) be fixed and let (u m ) m >o be the corresponding recursive 
sequence from ([9]). If n > 2 is an integer with gcd(n, D) = 1, ord(M, n) divides n(n). 
Moreover, with N n := n/ gcd(n, mgcd(M)) ; one has 

ord(M,n) = k(N u ) 

whenever N n > 1. In particular, this gives ord(M, n) = n(n) whenever n and mgcd(M) are 
coprime. 

In the remaining case, N n = 1, the matrix satisfies M = al mod n with a € (Z/nZ) x , so 
that ord(M, n) is the order of a modulo n. 

Proof. If M = ( " b d ) , the iteration formula (j!4j) implies that M m = 1 mod n if and only if 

u m a — Du m ^i = 1 , u m b = , u m c = , and u m d — Du m ^i = 1 mod n, 

so that also u m (a — d) = mod n. Consequently, n divides u m b, u m c and u m (a — d). This 
implies that u m is divisible by ged " ra;b) , gcd " niC) and ged(n n a _ d) , hence also by the least common 
multiple of these three numbers, which is the integer 

n n 
gcd(n, gcd(6, c, a — d)) gcd(n, mgcd(M)) 

Since N n \n, we now also have u m a — Du m -\ = 1 mod N n . When u m = mod N n , the 
recursion now gives u m+ i = Tu m — Du m -\ = —Du m -\ = l — u n a = 1 mod N n . Consequently, 
M m = 1 mod n is equivalent to u m = and w m +i = 1 mod N n . So, for N n > 1, one has 

ord(M,n) = K(N n ), 

which is the period of the sequence (u m )m>o modulo N n . Since k(N u ) clearly divides K(n), 
one finds ord(M, n) \ n{n). 

Finally, when N n = 1, one has n \ mgcd(M), which implies M = al mod n, where we have 
a 2 € (Z/nZ) x due to gcd(n,D) = 1. Since this also implies a € (Z/nZ) x , the last claim is 
clear. □ 
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Remark 1. Instead of the characteristic polynomial Pm, any other monic polynomial that 
annihilates M can be employed to derive a recursive sequence whose period is a multiple 
of the matrix order modulo n. For n = p a prime, the unique minimal polynomial Qm of 
M suggests itself to be chosen. For d = 2, Qm has smaller degree than Pm precisely when 
M = al, whence mgcd(M) = and Qm( z ) = z — a. Consequently, ord(M,p) is then always 
equal to the order of a modulo p. 

3. Orbit pretail structure of toral endomorphisms 

In this section, we look at the action of M € Mat(d, Z) on a lattice L n , with special emphasis 
on the structure of general endomorphisms. When M is not invertible, this manifests itself in 
the existence of non-trivial 'pretails' to periodic orbits, with rather characteristic properties. 
More precisely, given a periodic point y of M, a finite set of iterates (or suborbit) 

(15) O = {x, Mx, M 2 x, M n x = y} 

is called a pretail (of y) if y is the only periodic point of M in O. 

3.1. General structure. Let M and n be fixed, and define R = TLjnL. Let per(M) denote 
the set of periodic points on the lattice L n , under the action of M mod n. Due to the linear 
structure of M, per(M) is an M-invariant submodule of L n . It is the maximal submodule on 
which the restriction of M acts as an invertible map. The kernel kei(M k ) C L n denotes the set 
of points that are mapped to under M k . One has ker(M fc ) C ker(M fc+1 ) for all k > 0, and 
this chain stabilises, so that lJfc>o ker(M fc ) is another well-defined and M-invariant submodule 
of L n . This is then the maximal submodule on which the restriction of M acts as a nilpotent 
map. Note that per(M) n ker(M fc ) = {0} for all k > 0. 

Consider an arbitrary x € L n and its iteration under M. Since \L n \ = n d is finite, Dirichlet's 
pigeon hole principle implies that this orbit must return to one of its points. Consequently, 
every orbit is a cycle or turns into one after finitely many steps, i.e. it is eventually periodic. 
By elementary arguments, one then finds the following result. 

Fact 1. There are minimal integers m > and k > 1 such that M k+m = M m mod n. The 
number k is the least common multiple of all cycle lengths on L n , while m is the maximum 
of all pretail lengths. Clearly, per(M) = Fix(M fc ). □ 

The lattice L n = R d is a free R- module. The modules per(M) as well as Fix(M- ? ) and 
ker(M- J ) for j > 1 are submodules of it, with Fix(M*) nker(M J ) = {0} for alH > 1 and j > 0. 
Recalling some results on modules from |33|. Ch. Ill] now leads to the following consequences. 

Fact 2. Let m and k be the integers from Fact [TJ If m > 1, one has 

{0} C ker(M) C ker(M 2 ) C . . . C ker(M m ) C L n , 
while ker(M m+: ') = ker(M m ) for all j > 0. Moreover, one has 

L n = Fix(M fc ) ©ker(M m ), 
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Figure 1 . The directed graph for the action of M = ( f 4 ) on the lattice L 6 . 
The matrix has three fixed points and two 3-cycles (shown once each), while 
each periodic point has the same binary tree of height 2 as pretail tree. 

which is the direct sum of two M -invariant submodules. Hence, per(M) and ker(M m ) are 
finite projective R-modules. □ 

In general, the projective summands need not be free. As a simple example, let us consider 
L 6 with d = 1 and M = 2. Here, per(M) = {0, 2, 4} covers the fixed point and a 2-cycle, 
while ker(M) = {0,3}. Both are modules (and also principal ideals, hence generated by a 
single element) over Z/6Z, but do not have a basis, hence are not free. Nevertheless, one has 
Z/6Z = per(M) © ker(M). We will return to this question below. 

3.2. The pretail tree. Consider the equation M e x = y, with £ G N, for some arbitrary, but 
fixed y in L n . In general, this equation need not have any solution x € L n . On the other 
hand, when there is a solution x € L n , the set of all solutions is precisely x + ker(M ), which 
has cardinality |ker(M )|. If y is a periodic point, the first case can never occur, as there 
is then at least one predecessor of y. Due to the linearity of M, the structure of the set of 
pretails of a periodic point y must be the same for all y € per(M) (note that there is precisely 
one predecessor of y in the periodic orbit, which might be y itself, while all points of the 
pretail except y are from the complement of the periodic orbit). 

Consequently, we can study the pretail structure for y = 0. Let us thus combine all pretails 
of the fixed point into a (directed) graph, called the pretail graph from now on; see [UJ 
for general background on graph theory. A single pretail is called maximal when it is not 
contained in any longer one. By construction, there can be no cycle in the pretail graph, 
while y = plays a special role. Viewing each maximal pretail of as an 'ancestral line', we 
see that this approach defines a rooted tree with root 0. Note that an isomorphic tree also 
'sits' at every periodic point y. 

Corollary 1. Every periodic point of M on L n has a directed pretail graph that is isomorphic 
to that of the fixed point 0. Up to graph isomorphism, it thus suffices to analyse the latter. 
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Figure 2. The directed graph for the action of M = (5 1 £) on the lattice 
L 15 . The only fixed point of M is 0, while it has two 2-cycles and five 4-cycles 
(shown once each only). All pretail trees have the same height. 



By reversing the direction, it is a rooted tree with root 0. This tree is trivial if and only if M 
is invertible on L n . □ 

Two illustrative examples are shown in Figures [T] and [2j Each M defines a unique (rooted) 
pretail tree on a given lattice. If Uj denotes the number of nodes (or vertices) of this tree with 
graph distance i from the root, we have v = 1 and 

(16) | ker (M j ) \ = v + v 1 + ---+v j 

for all j > 0, where Vi = for all i larger than the maximal pretail length. Also, one has 

(17) Vj = |ker(M J ) \ker(M J ' _1 )| = | ker(Af j )| - | ker(M^ 1 )| 

for j > 1, where the second equality follows from the submodule property. Recall that 
terminal nodes of a rooted tree (excluding the root in the trivial tree) are called leaves. With 
this definition, the total number of leaves on L n is \L n \ ML n \. For i € N, define Wi to be the 
number of nodes with graph distance i from the root that fail to be leaves, and complete this 
with w = for the trivial tree and w = 1 otherwise. It is clear that this leads to 

(18) v i+1 = w {w i \\ B sr(M)\-6 i)0 ), 

via the number of solutions to Mx = and the special role of the root, and inductively to 

(19) | ker(M m )| = (w Q + w 1 + ■ ■ ■ + ^)| ker(M)| 
whenever u> = 1, with both relations being valid for all i > 0. 

Lemma 2. If M acts on L n , its uniquely defined pretail tree of the fixed point has height 
m > 0, and the following properties are equivalent. 

(i) All maximal pretails have the same length m; 

(ii) One has Vi = Wi ^ for all < i < m and Wi = for i > m; 

(iii) One has |ker(M m )| = |ker(M)| |ker(M i )| = |ker(M)| m for all < i < m, together 
with \ker(M m+ i)\ = |ker(M m )| for all j > 0. 
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In particular, m is the integer from Fact [TJ 

Proof. By Corollary [TJ the pretail tree is trivial (hence w = m = 0) if and only if M is 
invertible on L n . Since all claims are clear for this case, let us now assume that M is not 
invertible on L n . 

All maximal pretails have the same length if and only if all leaves of the pretail tree of 
have the same graph distance from the root 0. Clearly, the latter must be the height m of 
the tree. When M is not invertible on L n , the tree is not the trivial one, so m > 1. The 
equivalence of (i) and (ii) is then clear, since both conditions characterise the fact that all 
leaves have distance m from the root. 

The implication (ii) (iii) can be seen as follows. The first claim is trivial for i = 0, as 
ker(M°) = ker(l) = {0}. Assuming (ii), Eqns. ([IT ]) - ([15 ]) yield 

|ker(M i+1 )| = | ker(M i )| + Wi \ ker(M)| = | ker(AP)| + (|ker M l \ - |ker M !_1 1) | ker M| 

for 1 < i < m, which (inductively) reduces to the first condition of (iii), while the second is 
clear from the meaning of m. 

Conversely, the second condition of (iii) means w m +j = for all j > 0, while the first 
condition, together with Eqns. (|18p and (|19p . successively gives V{ = Wi for all < i < m. □ 

On the lattice L p r, when |ker(M)| = p, one can say more. 

Proposition 3. Consider the action of M on the lattice L p r. When |ker(M)| = p, one has 
|ker(AP)| = p mm (*> m ) for all i > 0, where m is the integer from Fact [T] for n = p r . This 
means Vi = p 1 (p — 1) for 1 < i < m, and all maximal pretails share the same length m. 

Proof. By assumption, M is not invertible, and the last claim is obvious from Lemma [2] in 
conjunction with Eq. (116[) . We thus need to prove the formula for the cardinality of ker(M*) 
for arbitrary i > 0. 

Since is a fixed point of M, we clearly have v Q = 1 and v-y = p — 1, together with 
the inequality < w-y < p — 1. If w-y = 0, we have m = 1 and we are done. Otherwise, 
ker(M) C ker(M 2 ), hence |ker(M 2 )| = p> for some j > 2, as the kernel is a subgroup of our 
lattice (which has cardinality p rd ). This forces = p — 1 and j = 2. More generally, when 
|ker(M*)| = p % for some 1 < i < m, one cannot have uii = 0, so that |ker(M* +1 )| = p % ^ 3 for 
some j > 1. Since now < w-i < p 4_1 (p — 1), the only possibility is j = 1 together with 
Wi = p l ~ l {p — 1). This argument can be repeated inductively until i = m is reached, with 
|ker(M m +J)| = |ker(M m )| for all j > 0. □ 

In general, the maximal pretails need not share the same length, which means that we still 
have to extend our point of view. 

Example 1. Consider the matrix M = (||) on L 8 , where it is nilpotent (mod 8) with 
nil-degree 4. Since card(ker(M)) = 4, Proposition [3] does not apply. The (directed) pretail 
graph spans the entire lattice and is shown in Figure [3] together with the loop at that marks 
this point as the root of the tree (which emerges from the figure by removing this loop and 
reversing all arrows). 
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(2,2) 



(2,0) 




(3,3) 











Figure 3. The pretail graph for Example with coordinates for the action 
of the matrix M on L s , where it is nilpotent with nil-degree 4. 



So far, we have looked at a single lattice L n . However, any given matrix M immediately 
defines a sequence of trees via L n with n € N. When d = 2, the result of [12, Thm. 2] implies 
the following result. 

Corollary 2. Let M,M' € Mat(2,Z) be two matrices with the same trace, determinant and 
mgcd. Then, they have the same sequence of pretail trees on the lattices L n . □ 

3.3. Decomposition on L p r. When the integ crs v ore coprime, one has Luv — ® 
wherefore the action on L n with n E N is completely determined by that on L p r, for all p r ||n. 
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In particular, the pretail orbit structure on an arbitrary L n can be derived from that on the 
sublattices associated with the factors in the prime factorisation of n. 

Define R r = r Ljp T 'L, which is a local ring, with unique maximal ideal (p) = pR r . The 
latter contains all zero divisors. By \33\ Thm. X.4.4], we then know that the two projective 
modules per(M) and kei(M k ) of Fact [2] are free, so each has a basis. Consequently, one 
knows that the linear map on L p r defined by M induces unique linear maps on Fix(M fc ( r )) 
and ker(M m ^ r - ) ), and M is conjugate to the direct sum of these maps, compare [TJ Prop. 4.3.28]. 
Each of the latter, in turn, admits a matrix representation with respect to any chosen basis 
of the corresponding submodule. Different choices of bases lead to conjugate matrices, by an 
application of [TJ Prop. 4.3.23]. 

Corollary 3. On L p r, M is similar to a block diagonal matrix (qb) over R r , where A is 
invertible and B is nilpotent, the latter of nil-degree n(B) say. The block matrices A and B 
are unique up to similarity. The direct sum from Fact [2] now reads 

L p r = Fix(M ord ( A ' p '')) ker(M n ( fi )), 

where the concrete form of the exponents k and m of Fact [2] follows from the block diagonal 
structure of M chosen. Here, Fix(M ord ( j4 ' pr )) ~ R d ' andker(M n ( B ^) ~ R d ~ d ' , where one has 
d' = rank (per (M)) < d. 

Furthermore, d' is independent of r. When comparing the above objects as modules over 
the ring R s for different s, one has 

rank 1 (per 1 (M)) = rank r (per r (M)) = d' and 
rank^ker^M" 1 ^))) = rank r (ker r (M m M)) = d - d' , 
where an index s at per, ker or rank refers to R s as the underlying ring. 

Proof. The diagonal block-matrix structure is clear from [l] Props. 4.3.28 and 4.3.23], while 
the isomorphism claim follows from [33, Cor. III. 4. 3]. 

For the last claim, observe that A and B can be viewed as integer matrices acting on R d ' 
and R d ~ d ' , respectively. Here, B s = modp r for some s G N and gcd(det(yl),p) = 1, because 
A is invertible mod p r and det(^4) must be a unit in R r . But this means that the reduction of 
A mod p is also invertible over R 1 = Z/pZ, while the reduction of B mod p is still nilpotent. 
Consequently, these reductions provide the blocks for the direct sum over L p , and the claim 
is obvious. □ 

Since two free modules of the same rank are isomorphic [331 Cor. III. 4. 3], we also have the 
following consequence. 

Corollary 4. One has the following isomorphisms of R\ -modules (as ¥ p -vector spaces), 
per r (M) jp per r . (M) ~ per x (M) and ker r (M m ^ ) /p ker r (M m ^ ) ~ kerf (M m W ) . 
This implies 

| per r (M) | = p rd ' = \pex 1 (M)\ r and \ ker r (M m ^ ) | = p r ^- d '^ = \ k eri (M m « ) | r 
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for the cardinalities of the finite modules. □ 

At this point, it is reasonable to link the properties of M on L p r to its minimal polynomial 
over W p . 

Lemma 3. If M is similar mod p to the block diagonal matrix of Corollary O its minimal 
polynomial over ¥ p is fJ> M (x) = x s f(x), where f is a monic polynomial of order k over ¥ p 
with /(0) 7^ 0. When M is invertible, one has s = and k = gcd{£ € N | M l = 1 modp}. 
When M is nilpotent, f = 1 and s = gcd{t £ N | M* = mod j>}. In all remaining cases, s 
and k are the smallest positive integers such that B s = and A k = 1 mod p. 

Proof. Recall from [344 Def. 3.3.2] that the order of a polynomial / € ¥ p [x] with /(0) ^ 0, 
denoted by ord(/,p), is the smallest positive integer £ such that f(x)\(ar — 1). When M is 
invertible and k as claimed, the polynomial x k — 1 annihilates M. Since f^ M (0) ^ in our case, 
we have [i M = f with f{x)\{x k — 1), so that ox&(f,p)\k by [Ml Lemma 3.3.6]. By construction, 
k is also the minimal positive integer such that x k — l annihilates M, hence k = ovd(f,p). 

When M is nilpotent, the claim is obvious, because is then the only possible root of the 
minimal polynomial over ¥ p , as all other elements of the splitting field of / are units. 

In all remaining cases, M is similar to A © B with A invertible and B nilpotent, by 
Corollary [3l We thus know that n M (x)\x s (x k — 1) with s and k as claimed, since the latter 
annihilates both A and B. Observe that B s (B k — 1) = mod p means B k+S = B s mod p. 
Since B is nilpotent, its powers cannot return to a non-zero matrix, hence B s = mod p. 
Similarly, A s+k = A s mod p is equivalent with A k = 1 mod p, as A is invertible. This shows 
that we must indeed have n M (x) = x s f{x) with ord(/,p) = k. □ 

3.4. Classification on L p . When we consider n = p, we can go one step further, because 
F p is a field and one can classify nilpotent matrices via their Jordan normal form. This 
follows from the observation that is the only possible eigenvalue. Recall that an elementary 
shift matrix is an upper triangular matrix with entries 1 on the upper super-diagonal and 
everywhere else (this includes the 0- matrix in one dimension). An elementary shift matrix 
is nilpotent, with nil-degree equal to its dimension. The following result is now a standard 
consequence of the Jordan normal form over fields [271 [33] . 

Fact 3. The nilpotent matrices in Mat(d, F p ) are conjugate to block- diagonal matrices, where 
each block is an elementary shift matrix. □ 

Some of this structure survives also for general n. For instance, the 0-matrix in dimension 
d > 1 leads to the regular (n d — l)-star as its pretail tree on L n . When d>2, the d-dimensional 
elementary shift matrix, on L n , results in a semi-regular tree with w = 1, w l = n — 1 and 
Wi = n for 2 < i < d — 1, while Wj = for all j > d. These trees have the property that all 
maximal pretails share the same length, which is the nil-degree of the matrix. One can now 
go through all possible block-diagonal combinations of such elementary shift matrices. This 
is a combinatorial problem and gives the possible pretail trees over ¥ p . 
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As already suggested by Proposition [3l the structure of ker(M) plays an important role for 
the structure of the pretail tree. Together with the linearity of M, it constrains the class of 
trees that are isomorphic to the pretail tree of some integer matrix. A more detailed analysis 
is contained in |36j . 

4. Symmetry and reversibility 

Reversibility is an important concept in dynamics, compare [39] and references therein for 
background, and [20] for an early study in continuous dynamics. Here, we focus on discrete 
dynamics, as induced by toral auto- and endomorphisms. 

A matrix M is called reversible, within a given or specified matrix group Q, if it is conjugate 
to its inverse within Q. Clearly, this is only of interest when M 2 ^ 1. To put this into 
perspective, one usually defines 

S(M) = {GeG\ GMG- 1 = M} and 1Z{M) = {G € Q \ GMG' 1 = M ±l } 

as the symmetry and reversing symmetry groups of M; see and references therein for 
background and [9j [10] for examples in our present context. In particular, one always has 
1Z(M) = S(M) when M 2 = 1 or when M is not reversible, while 1Z(M) is an extension of 
S{M) of index 2 otherwise. 

Note that a nilpotent matrix M (or a matrix with nilpotent summand, as in Corollary [3]) 
cannot be reversible in this sense. However, they can still possess interesting and revealing 
symmetry groups, although it is more natural to look at the ring of matrices that commute 
with M in this case. 

Example 2. Reconsider the matrix M = (f |) from Example Q3 and its action on L 8 . 
Clearly, M commutes with every element of the ring Z/8Z [M], which contains 64 elements. 
This follows from the existence of a cyclic vector, but can also be checked by a simple direct 
calculation. Consequently, the symmetry group (in our above sense) is the intersection of this 
ring with GL(8, Z/8Z), which results in 

S{M) = ((if) ,3-1,5-1) ~ C 8 x C 2 x C 2 , 

which is an Abelian group of order 32. The matrices in S(M) have either determinant 1 or 
5, with {,4 G S(M) | det(^) = 1} ~ C 4 x C 2 x C 2 . 

One can now study the action of S(M) on the pretail graph of Figure O which actually 
explains all its symmetries. 

In what follows, we derive certain general properties, where we focus on the reversing 
symmetry group, with invertible matrices M in mind. 

4.1. Reversibility of SL(2, Z)-matrices mod n. Recall the matrix mgcd from Eq. (|13[) . 
which is a conjugation invariant. It can be used to solve the reversibility at hand as follows. 

Theorem 1. Let M € SL(2,Z) and n € N be arbitrary. Then, the reduction of M mod n is 
conjugate to its inverse within the group GL(2,Z/nZ). The action mod 1 of any M E SL(2,Z) 
on L n is thus reversible for all n G N. 
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Moreover, if M G SL(2, Z) has mgcd(M) = r ^ 0, its reduction mod n, for every n G N, 
possesses an involutory reversor. 

Proof. When M € SL(2,Z), also its inverse is in SL(2,Z), and M and M~ l share the same 
determinant and trace. Moreover, they also have the same mgcd, so that the first claim 
follows from \12\ Thm. 2] (or from Lemma [1]). This immediately implies, for all n G N, 
the reversibility of the action mod n of M on the lattice L n , so that the statement on the 
equivalent action of M mod 1 on L n is clear. 

Now, let M = [ a c b d ) G SL(2,Z), so that M _1 = ( _? c and M and M" 1 share the same 
determinant (1), trace (a + d) and mgcd (r). Assume r 7^ 0, let n > 2 be fixed and consider 
the matrices mod n. Recall the normal forms 

N(M) = I T | and NiM' 1 ) = 
\r d I 

as defined in the proof of [12} Prop. 6], and note that they are not inverses of each other. 
However, by [12 Prop. 5], there is some matrix P n G GL(2,Z/nZ) with M = P n N(M)p- 1 , 
hence we also have M _1 = P n (N{M)) P~ l . Ob serve next that 




(N(M)y 1 = ( d J ) = c[l I ) C" 1 = CNiM-^C- 1 , 



be 





where C = (0-1) is an involution. On the other hand, N(M) and N(M 1 ) satisfy the 
assumptions of \12\ Prop. 6], so that 

N(M~ 1 ) = AN(M)A^ with A 

where we globally have A = ( -1 ) whenever d = a in the original matrix M. Together with 
the previous observation, this implies (iV(M)) 1 = (CA)N(M)(CA)~ 1 where 

CA = 

is an involution. Putting everything together, we have 

M" 1 = {P n (CA)p- 1 )M(P n (CA)p- 1 y\ 
which is the claimed conjugacy by an involution (which depends on n in general). □ 

Note that the matrix M in Theorem [T] need not be reversible in GL(2,Z), as the example 
M = (f f 6 ) from [9, Ex. 2] shows. Nevertheless, for any M G SL(2,Z) with mgcd(M) / 
and n > 2, the (finite) reversing symmetry group of M within GL(2,Z/nZ) is always of the 
form 7Z(M) = S(M) x C2, with C2 being generated by the involutory reversor. The structure 
of S{M) remains to be determined. 

In the formulation of Theorem [TJ we have focused on matrices M G SL(2,Z) because the 
condition tr(M) = tr(M _1 ) for a matrix M with det(M) = —1 forces tr(M) = 0, which means 
that M is itself an involution (and thus trivially reversible in GL(2,Z)). More interesting 
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(beyond Theorem [T]) is the question which matrices M £ Mat(2,Z), when considered mod n 
for some n € N, are reversible in GL(2, Z/nZ). Let us begin with n = p being a prime, where 
Z/pZ ~ F p is the finite field with p elements. 

4.2. Reversibility in GL(2,F p ). Let us consider the symmetry and reversing symmetry 
group of an element of GL(2,F p ) with p prime, the latter being a group of order 

|GL(2,F P )| = (p 2 -l)(p 2 - p ) = p(p-l) 2 (p + l), 

compare Eq. (|3j). For our further discussion, it is better to distinguish p = 2 from the odd 
primes. For convenience, we summarise the findings also in Table 14.21 

Example 3. For p = 2, one has GL(2,F 2 ) = SL(2,F 2 ) ~ D 3 , the latter denoting the 
dihedral group of order 6. There are now three conjugacy classes to consider, which may be 
represented by the matrices 1 , the involution R = ( 5 o ) > ana - ^ ne ma trix M = ( \ q ) °f order 
3. The corresponding cycle structure on L 2 is encapsulated in the generating polynomials 
Z 2 {t). They read 

(1-i) 4 , (l-t) 2 (l-t 2 ) and (l-t)(l-t 3 ), 

respectively, and apply to entire conjugacy classes of matrices. 

For the (reversing) symmetry groups, one clearly has 7^(1) = 5(1) = GL(2,F 2 ), while 
1Z{R) = S(R) = (R) ~ C 2 . The only nontrivial reversing symmetry group occurs in the third 
case, where S{M) = (M) ~ C 3 . Since RMR = M 2 = M" 1 , one has K(M) = GL(2,F 2 ) ~ 
C3 x C 2 . So, all elements of GL(2,F 2 ) are reversible, though only M and M 2 are nontrivial 
in this respect. 

For p an odd prime, one can use the normal forms for GL(2, F p ), see [331 Ch. XVIII. 12], to 
formulate the results; compare Table [4T2l We summarise the reversibility and orbit structure 
here, but omit proofs whenever they emerge from straight-forward calculations. 

I. The first type of conjugacy class is represented by matrices M = at with a G F^ ~ C p -\. 
The order of M coincides with the order of a mod p, ord(a,p), which divides p — 1. One clearly 
has K(M) = S(M) = GL(2,F P ) in this case, either because a 2 = 1 (so that M = M' 1 ) or 
because a 2 ^ 1 (so that no reversors are possible). The corresponding orbit structure on 
L p comprises one fixed point (x = 0) together with f d (~p) orbits of length ord(a,p). The 
non-trivial orbits starting from some i / must all be of this form, as x gets multiplied by 
a under the action of M and returns to itself precisely when a k = 1, which first happens for 
k = ord(a,p). 

II. The next type of conjugacy class is represented by matrices M = (g with a £ W* . 
Its symmetry group is given by 

which is Abelian. As generators of the cyclic groups, one can choose (J J), which has order 
p in GL(2, F p ), and 7I, with 7 a generating element of F^ . The reversible cases are precisely 
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the ones with a 2 = 1 in F p , hence with det(M) = 1. Here, R = diag(l,— 1) is a possible 
choice for the (involutory) reversor, so that TZ(M) = S(M) x (R) ~ (C p x C p _i) x Ci. 
A matrix M of type II (in its normal form as in Table I4.2p satisfies 

whence a point (x,0) with x ^ is fixed by M k if and only if k = ord(a,p), and a point 
(x,y) with xy 7^ if and only if p\k and ord(a,p)|fc. Since ord(a,p)|(p — 1), one has 
lcm (p, ord(a,p)) = 1, wherefore this gives OT ^^ p ^ orbits of length p — 1 and p—fe^ = ord^p) 
orbits of length pord(a,p) in total. 

III. The third type of conjugacy class is represented by M = diag(a, b) with a, & £ F x 
and a ^ b. This results in S(M) = {diag(a, /3) | a, f3 G F*} ~ Cp_ 1 . The condition for 
reversibility leads either to a 2 = b 2 = 1, hence to 6 = —a, or to ab = 1. In the former 
case, M itself is an involution, so that 1Z(M) = S(M) is once again the trivial case, while 
det(M) = ab = 1 leads to genuine reversibility, with involutory reversor R = (5 o) an< ^ hence 
to K{M) = S{M) x C 2 . 

For a type III matrix, one has M k (x,y) t = (a fc x,6 fc y)*, so each of the p — 1 non-zero 
points (x, 0)* is fixed by M ord ( a ' p ) ; analogously, each of the p — 1 non-zero points (0, y)* 
is fixed by M ord ( b ' p ). The remaining points that are non-zero in both coordinates have 
period lcm(ord(a,p), ord(b,p)). In summary, this gives one fixed point, 3^7^^ orbits of 

length ord(a,p), ^f^y orbits of length ord(&,p), and lcm{o J(~p^ mdM) orbits of length 
lcm(ord(a,p), ord(6,p)). 

IV. Finally, the last type of conjugacy class can be represented by companion matrices 
of the form ( 9 ~® ) with the condition that the characteristic polynomial z 2 — Tz + D is 
irreducible over F p . The determinant and the trace satisfy D = rjrj' and T = 77 + 77', where r\ 
and rj' are not in ¥ p , but distinct elements of the splitting field of the polynomial, which can 
be identified with F p 2. One consequence is that 1 + D ± T = (1 ± r/)(l ± 77') 7^ 0. 

The symmetry group is 5(M) = {al + 7M | a, 7 G F p , not both 0}, which is an Abelian 
group with p 2 — 1 elements. The order follows from the observation that det(al + 7M) = 
(a + 777) (a + 777') vanishes only for a = 7 = in this case. In fact, one has S{M) ~ C p 2_ 1 , 
as any matrix ^ j G GL(2,F p ) with 77 G F p 2 \ F p has order p 2 — 1 or possesses a root in 
GL(2,F p ) of that order. This relies on the facts that we can always write 77 = A" 1 , where A is 
a generating element of F* 2 ~ C p 2_i, and that AA' and A + A' are in F p . This is a special case 
of Fact [9] below and of a statement on the existence of roots in GL(d, Z); see Lemma [6] below. 

The condition for reversibility, in view of the above restriction on D and T, can only be 
satisfied when D = 1, in which case R = (5 ) turns out to be an involutory reversor, so that 
again 1Z(M) = S(M) x C2 in this case. 

Matrices with irreducible characteristic polynomial \m P r °d uce orbits of one length r only, 
where r is the smallest integer such that Xm( z )\( zT ~ 1)j or ' equivalently, the order of its roots 
in the extension field F ,. 
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Table 1. Summary of conjugacy structure for GL(2,F p ) via normal forms. 
Note that class III is absent for p = 2. The second possibility for 1Z(M) always 
applies when det(M) = 1. Only non-trivial orbits are counted. 



class 


I 


II 


III 


IV 


normal form 
of matrix class 


al 

„ r- I?X 

a t ffp 


( a 1\ 
\ a 1 

a t JTp 


U b) 

r, -J- h CZ W X 
a ^ t ifp 


1 —D \ 

z 2 — '/ 2 -)- irred. 


min. polynomial 


(z — a) 


(z - a) 2 


(z — a)(z — 6) 


z 2 -Tz + D 


size of class 


1 


p 2 - 1 




2 

P — P 


no. of classes 


p — 1 


p — 1 


i(p_l)(p_2) 


\p{p-l) 


S(M) 


GL(2,F P ) 


Cp x Cp-i 


Cp-i x Cp-i 


Cp 2 -i 


1Z(M) 


S{M) 


5(M) or 
«S(M) x C 2 


5(M) or 
S(M) x C 2 


5(M) or 
5(M) x C 2 


orbit length 


ord(a,p) 


see text 


see text 


ord(x M ,p) 


orbit count 


p 2 -i 


see text 


see text 


p 2 -i 


ord(a,p) 


ord (X M 'P) 



Putting these little exercises together gives the following result. 

Theorem 2. A matrix M £ GL(2, F p ) is reversible within this group if and only if M 2 = 1 
or det(M) = 1. Whenever M 2 = 1, one has H(M) = S(M). If det(M) = 1 with M 2 ^ t, 
there exists an involutory reversor, and one has 1Z(M) = S(M) x C 2 . □ 

Remark 2. Since ¥ p is a field, we can use the following dichotomy to understand the structure 
of S(M), independently of the chosen normal forms. A matrix M € GL(2,F p ) is either 
a multiple of the identity (which then commutes with every element of Mat(2,F p )) or it 
possesses a cyclic vector (meaning an element v G F 2 such that v and Mv form a basis of F 2 ). 
In the latter case, M commutes precisely with the matrices of the ring F p [M], and we have 
S(M) = ¥ p [M] x = F p [M] nGL(2,F p ). This systematic approach provides an alternative (but 
equivalent) parametrisation of the above results for the normal forms. 

The question for reversibility in GL(2,Z/nZ) with general n is more complicated. The 
matrix M = (J ~^ 4 ) is reversible over Z/3Z (where it is an example of type IV), but fails to 
be reversible over Z/9Z, as one can check by a direct computation. Here, zero divisors show 
up via non-zero matrices A with AM = M _1 A, but all of them satisfy det(A) = mod 9. In 
fact, one always has A(Lg) C L% here. 

In general, the relation AM A' 1 = M' 1 with A, M € GL(2,Z/nZ) implies MAM = A 
and hence det(M) 2 = 1, because det(A) G (Z/nZ) x . Over F p , this gives det(M) = ±1, with 
reversibility precisely for det(M) = 1 according to Theorem [2j In general, one has further 
solutions of the congruence m 2 = 1 mod n, such as m = 3 for n = 8 or m = 4 for n = 15. 
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In any such case, M = J") is a matrix with M 2 = —ml. When m ^ — 1 mod n, M is 
of order 4 in GL(2,Z/nZ). It is easy to check that RMR = M' 1 = ( _° m I) in GL(2,Z/nZ), 
with the involution i? = (9 g). This establishes reversibility with 7Z(M) = S(M) xi C2. 

4.3. Some extensions to higher dimensions. In principle, a similar reasoning, based on 
a normal form approach, can be applied to arbitrary dimensions. Over the finite field ¥ p , 
normal forms are given by the rational canonical form and the elementary divisor normal 
form ('First' and 'Second natural normal form' in the terminology of [241 §6]), which are 
block diagonal matrices with companion matrices on the diagonal. 

The advantage of dealing with companion matrices is that one can employ the theory of 
linear recursions: there is a one-to-one correspondence between the cycle lengths modulo 
n € N of a certain initial condition u = (no, • • • ,Ud-i) under the recursion induced by the 
polynomial /, and the period of the corresponding point iz* under the matrix iteration of Cy; 
compare the final remark in |48| , and Section 12.41 

Working with a block diagonal matrix of this shape, the analysis can be done block-wise; 
in particular, the symmetry groups are the direct product of the symmetry groups of the 
component matrices on the diagonal, augmented by all additional symmetries that emerge 
from equal blocks, which can be permuted. 

Determining the period lengths associated with irreducible polynomials amounts to finding 
their orders in the sense of |34t Def. 3.3.2]. The periods and their multiplicities arising from 
the powers of irreducible polynomials that show up in the factorisation of the invariant factors 
(the elementary divisors) are then given by [481 Thm. 4]. 

Extending the analysis to matrices over the local rings Z/p r Z is more difficult. In general, 
it seems hard to write down an exhaustive system of normal forms for the similarity classes, 
and to decide whether given matrices are similar. However, a solution for a large subclass 
of square matrices over the p-adic integers Z p and the residue class rings Z/p r Z is presented 
in [TTj. For a polynomial / £ Z p [x] whose reduction modulo p has no multiple factors, a 
complete system of d x d matrix representatives X with respect to similarity that satisfy 
f(X) = mod p r is given by all direct sums of companion matrices which are in agreement 
with the factorisation of / mod p. For instance, if the reduction of the common characteristic 
polynomial modulo p of two matrices does not have any quadratic factors, the matrices are 
conjugate mod p T if and only if they are conjugate mod p \17\ Thm. 3 and Corollary]. 

An exhaustive treatment of conjugacy classes of 3 x 3 matrices over an arbitrary local 
principal ideal ring can be found in [3]. 

Remark 3. In [4J, it is pointed out that 2x2 matrices over a local ring can be decomposed 
into a scalar and a cyclic part. Over Z/p r Z, this decomposition reads 

M = dl + p £ C, 

where = gcd(mgcd(M), p r ) = pM m § cd ( M )) with v p denoting the standard p-adic valuation, 
unique d € {5_w=o a jP'' I P \ a j} an d cyclic C € Mat (2, Z/p r_f Z), which is unique up to 
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similarity. Moreover, C can be chosen as a companion matrix with the appropriate trace and 
determinant. 

Since dl and C commute, powers of M can be expanded via the binomial theorem. Using 
that the binomials satisfy gcd /^ k \ | (^) , the period per(x,p r ) of all x £ L p r is bounded by 

pei(x,p r ) < oid(d,p r ) ■ p r ~ , 

provided that 1 < I < r. Let ELf : Z/p r Z — > Z/p'Z denote the canonical projection, and let 
Sj(A) be the symmetry group of an integer matrix A, viewed as a matrix over TLjpPTL. Then, 
for p ^ 2 and £ > 1, one obtains S r (M) = ILJ 1 (Se(C)) from the symmetry equations. 

4.4. Reversibility mod n. Let M be a general integer matrix, with determinant D. 

Fact 4. 7/ M £ Mat(d, Z) is reversible mod n, one has D 2 = 1 mod n. Moreover, reversibility 
for infinitely many n implies D = 1 or D = — 1. 

Proof. The reversibility equation yields detM = detM -1 , hence D 2 = 1 mod n. If D 2 — 1 
has infinitely many divisors, one has D 2 = 1, hence D = lorZ) = — 1. □ 

Before we continue with some general result, let us pause to see what Fact H] specifically 
implies for d = 2. 

Fact 5. If M £ Mat(2,Z) with D = — 1 mod n is reversible mod n, one has 2tr(M) = 
mod n. In particular, tr(M) = mod n holds whenever n is odd. 

Proof. The trace is a conjugacy invariant, so reversibility mod n implies tr(M) = tr(M _1 ) 
mod n. The inversion formula for 2x2 matrices yields tr(M _1 ) = ^Ml = — tr(M) mod n, 
and thus 2 tr(M) = mod n. □ 

Fact 6. Consider M £ Mat(2,Z) with D = —1 mod n. Then, M is an involution mod n if 
and only if tr(M) = mod n. 

Proof. Let M = (" With D = —1, the inversion formula for M shows that M = M~ 1 is 
equivalent to d = -a. Thus, M 2 = 1 if and only if tr(M) = 0. □ 

The previous two facts imply 

Corollary 5. Let M £ Mat(2,Z) be reversible mod n > 2 with D = — 1 mod n. Then, 
M 2 = 1 mod n for n odd, and M 2 = 1 mod n/2 for n even. □ 

Let us continue with the general arguments and formulate a necessary condition for local 
reversibility. 

Lemma 4. Let p / 2 be a prime. If M € Mat(d, Z) is reversible mod p r , one has D = ±1 
mod p r . If d = 2, M is reversible mod p r if and only if D = 1 or M 2 = 1 mod p r . 

If M £ Mat(d, Z) is reversible mod 2 r , then D = ±1 mod 2 r_1 . When d = 2 and M is 
reversible with D = —1 mod 2 r ~~ 1 , one has M 2 = 1 mod 2 r ~ 2 . 
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Proof. For p ^ 2, Fact [H implies D 2 = 1 mod p r . Since p cannot divide both D — l and D + 1, 
one has p r \(D — 1) or p r \(D + 1), which gives the first claim. When 2 r \(D — 1)(D + 1), 2 
divides one of the factors and 2 r_1 the other one, so D = 1 or D = — 1 mod 2 r . If D = — 1 
mod 2 r ~\ Fact Ogives 2 tr(M) = mod 2 r " 1 and thus M 2 = 1 mod 2 r " 2 by Fact El □ 

One immediate consequence for d = 2 is the following. 

Corollary 6. If M E GL(2,Z) wii/i D = —1 is reversible for infinitely many ra E N, one has 
M 2 = t. □ 

Fact 7. Lei ^4 5e an integer matrix whose determinant is coprime with n E N. The reduction 
of the inverse of A over Z/nZ, taken modulo k\n, is then the inverse of A over "L/kZ,. □ 



Lemma 5. Let n = p 7 ^ 1 . . .p r s s be the prime decomposition of n E N. Then, two matrices 
M,M' E Mat(d, Z) are conjugate mod n if and only if they are conjugate mod p r ^ for all 
l<i<s. 

Proof. M ~ M' mod n means M' = AMA^ 1 for some A E GL(n, Z), which implies conjugacy 
mod k for all k\n. 

For the converse, let Ai E GL(d, Z/p[*Z) denote the conjugating matrix mod p[\ The Chi- 
nese remainder theorem, applied to each component of the matrices Ai and A' 1 , respectively, 
gives matrices A and B that reduce to A{ and A^ modulo p\\ respectively. By construction, 
AB = 1 mod pY for all i, hence also AB = 1 mod n and thus B = A^ 1 in GL(d, Z/nZ). □ 

Proposition 4. VFii/i n as in Lemma® a matrix M E Mat(d, Z) is reversible mod n if and 
only if M is reversible mod p^ for all 1 < i < s. 

Proof. The claim is a statement about the conjugacy of M and M _1 in the group GL(d, Z/nZ), 
which is thus a consequence of Lemma [5l We just have to add that, by Fact [71 the inverse of 
M mod n reduces to the inverse mod p^ 1 , so MR = RM^ 1 mod p r ^ for all i. □ 

Corollary 7. Consider a matrix M E Mat(2,Z) mi/i D = det(M) and let n = p r ^p r 2 ■ ■ -Ps s - 
When n is not divisible by A, M is reversible mod n if and only if, for each 1 < i < s, D = 1 
or M 2 = 1 mod p^ . When n = 2 ri p2 2 • • - Pi" with r 1 > 2, M is reversible mod n if and only 
if it is reversible mod 2 ri and, for all i > 1, D = 1 or M 2 = 1 mod p^ . 

Proof. According to Lemma [5l the matrix M is reversible mod n if and only if it is reversible 
mod for all 1 < i < s. By Lemma HI this is equivalent with D = 1 or M 2 = 1 mod for 
alii with A\p\\ □ 

Remark 4. To see that reversibility mod p for all primes p which divide n is not sufficient for 
reversibility mod n, one can consider a locally reversible matrix M with det M/l: according 
to Fact(H only finitely many n exist such that M is reversible mod n, so for each prime p there 
must be a maximum r for which M is reversible mod p r . Recalling an example from above, 
M = ( 5 ~o 4 ) is reversible mod 3 but not mod 9 as can be verified by explicit calculation. It 
is an involution mod 5, hence also reversible mod 15, but not mod 45. 
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Reversibility can be viewed as a structural property that reflects additional 'regularity' in 
the dynamics, in the sense that it typically reduces the spread in the period distribution. For 
2 x 2-matrices, the normal form approach shows that reversibility implies the existence of 
only one non-trivial period length on L p ; compare our comments in Section [4.51 

4.5. Matrix order and symmetries over F p . Let us now discuss the order of a matrix 
M G Gh(d, F p ), with p a prime, in conjunction with the existence of roots of M in that group. 
We begin by recalling the following result from |34|, Thm. 2.14, Cor. 2.15 and Cor. 2.16]. 

Fact 8. If f is an irreducible polynomial of degree d over ¥ p , its splitting field is isomorphic 
with F p d. There, it has the d distinct roots a,a p , . . . ,a pd 1 that are conjugates and share the 
same order in (F p d) x . 

In particular, two irreducible polynomials over ¥ p of the same degree have isomorphic 
splitting fields. □ 

From now on, we will identify isomorphic fields with each other. In particular, we write 
¥ p d for the splitting field of an irreducible polynomial of degree d over F p . 

Next, let K be an arbitrary finite field, consider an irreducible, monic polynomial / G K[x] 
of degree d, and let L be the splitting field of /. When A 1; A 2 , • • • , X d are the roots of / in L, 
one has the well-known factorisation 
d 

(20) f{x) = Hix-Xj) = x d -e l (X 1 ,...,X d ) + ... + (-l) d e d (X 1 ,...,X d ), 

3=1 

where the denote the elementary symmetric polynomials, 

e i( x ii •••> x d) = x i + x 2 + •••+ x d , ... , e d (x ll . . . , x d ) = x 1 ■ x 2 ■ ■ ■ ■ ■ x d . 

The elementary symmetric polynomials, when evaluated at the roots of /, are fixed under all 
Galois automorphisms of the field extension L/K, so that the following property is clear. 

Fact 9. An irreducible, monic polynomial f G K[x] satisfies (|20p over its splitting field L. 
In particular, the elementary symmetric polynomials e 1 ,...,e d , evaluated at the d roots of f 
in L, are elements of K . □ 

Let M be a d x d integer matrix with irreducible characteristic polynomial Xm over ^p- 
Let a be a root of Xm ^ n ^p d an< ^ X a generating element of the unit group (F p d) x . Clearly, 
there is an n G N with a = X n . By Fact [SJ one has W p (a) = ¥ p d = F p (A), where the degree of 
the extension field over ¥ p equals d. Consequently, the minimal polynomial of A over F p is an 
irreducible monic polynomial of degree d over F p , and the conjugates of a are powers of the 
conjugates of A. Let a 1 , . . . ,a d and X 1 , . . . ,X d denote the respective collections of conjugates. 
Thus, over ¥ p d, one has the matrix conjugacy 

M ~ diag(a 1 ,...,a d ) = diag(A ls . . . , X d ) n ~ C(f) n , 

with f(x) G ¥ p [x] as in (|20|) and C(f) denoting the companion matrix of /. Here, it was 
exploited that a d x d matrix whose characteristic polynomial / has d distinct roots is always 
similar to the companion matrix of /. Note that C(f) G GL(d,F p ) by Fact [9l 
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Now, M and C(f) are matrices over F p that are conjugate over ¥ p d, so (by a standard 
result in algebra, see [TJ Thm. 5.3.15]) they are also conjugate over F„, which means that we 
have the relation 

(21) M = A~ x C(f) n A = (A- l C(f)A) n =: W n 

with some A G GL(d,F p ). By similarity, ord(IF) = ord(C(/)) = ord(diag(A 1 , . . . , A d )) = 
p d — 1. This gives the following result. 

Lemma 6. A matrix M G GL(d, F p ) with irreducible characteristic polynomial has either the 
maximally possible order p d — 1, or admits an n-th root W £ GL(d,¥ p ) as in (|2ip . Here, n 
can be chosen as n = ^jpjy, so that the root has order p d — 1. □ 

Fact 10. Let A be a matrix over ¥ p with minimal polynomial of degree d. Then, the ring 

¥ P [A] = {^i + .^ + e^- 1 \i 3 e¥ p } 
has precisely p d elements, which correspond to the different d-tuples (£i, . . . , £d)- 

Proof. Two distinct d-tuples producing the same matrix would give rise to a non-trivial linear 
combination that vanishes, involving powers of A of degree d — 1 at most, which contradicts 
the minimal polynomial having degree d. □ 

Lemma 7. Let W,M € GL(d,¥ p ) satisfy W n = M and ord(W r ) =p d -1. Then, ¥ p [M] = 
¥ P [W] and 

F P [M] X = F p [M]\{0} = (W) ~ C pd _ t , 
where (W) denotes the cyclic group generated by W. 

Proof Clearly, ¥ P [M] = ¥ p [W n ] C ¥ p [W], while Fact M implies \¥ p [M]\ = \¥ p [W]\ = p d , 
whence we have equality. Further, 

(W) C ¥ P [W] X C ¥ P [W}\{0} = F p [M]\{0}, 

and again, comparing cardinalities, one finds |(W)| = p d - 1 = \¥ p [M] \ {0}|, from which the 
claim follows. □ 

Let us summarise and extend the above arguments as follows. 

Corollary 8. A dxd integer matrix M with irreducible characteristic polynomial over the 
field ¥ p has a primitive root W £ GL(ci, F p ) with ord(W) = p d — 1. Moreover, one then has 
¥ p [M] x = ¥ p [M] \ {0} = (W) ~ C pd _ x . In particular, S{M) ~ C pd _ 1 in this case. 

More generally, we have S{M) = F P [M] X whenever the minimal polynomial has degree d. 

Proof. Since we work over the field F p , the irreducibility of the characteristic polynomial of 
M means that the minimal polynomial agrees with the characteristic polynomial and has 
thus maximal degree d. This situation is equivalent with M being cyclic [27, Thm. III. 2]. By 
Thm. 17 of |27| and the Corollary following it, we know that any matrix which commutes 
with M is a polynomial in M, so that S(M) = F p [M] x is clear. 
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The claim for matrices M with an irreducible characteristic polynomial follows by Lemmas[6] 
and CD □ 

When a matrix M G Mat(cZ, ¥ p ) fails to be cyclic, there are always commuting matrices 
that are not elements of F p [M], see Thm. 19 of |27j and the following Corollary. In such a 
case, S{M) is a true group extension of F P [M] X . The situation is thus particularly simple for 
matrices M G Mat(2,F p ): Either they are of the form M = at (then with S (M) = GL(2,F P )), 
or they are cyclic (then with S(M) = F p [M] x ). 



Appendix: Two classic examples 

If one reads through the literature, two matrices are omnipresent as examples, the Arnold 
and the Fibonacci cat map. Still, several aspects of them are unclear or conjectural, despite 
the effort of many. Let us sum up some aspects, with focus on properties in line with our 
above reasoning. 

A.l. Arnold's cat map. Here, we collect some results for the matrix Ma = ( \ \ ) G SL(2, Z) 
in an informal manner. This case was studied in [37\ |2T1 [25] and appeared in many other 
articles as main example. It was introduced in [3], Example 1.16] as a paradigm of (discrete) 
hyperbolic dynamics. 

The integer matrix Ma is reversible within the group GL(2,Z), with a reversor of order 
4, but none of order 2. One has S(M\) — C2 x Coo, where C2 = {±1} and the infinite 
cyclic group is generated by the unique square root of Ma in GL(2,Z) (see below), while 
TZ(Ma) = <S (A/a) x C4; see [9] for more. In particular, Ma inherits local reversibility in 
GL(2,Z/pZ) for all primes p from its 'global' reversibility within GL(2,Z). 

It was shown in [25J that Ma, except for the trivial fixed point 0, has orbits of only one 
period length on each prime lattice L p . In view of the normal forms, this is clear whenever 
the characteristic polynomial is irreducible. However, a matrix of type III from Table 14.21 
has reducible characteristic polynomial and occurs for primes with (|) = —1. Here, different 
orbit lengths would still be possible in general, but reversibility forces the two roots to be 
multiplicative inverses of one another and thus to have the same order modulo p. 

The iteration numbers are p m = f2m, where the ft are the Fibonacci numbers, defined by 
the recursion fk+i = ft + fk-i for A; G N with initial conditions /o = and f\ = 1. Since 
mgcd(M0 = 1; Proposition [2] implies 

ord(M A ,n) = « A (n) = periodj (/ 2 m)m>o mod n}, 

where the periods for prime powers (with r G N) are given by 

K A (2 r ) = 3 • 2 max ^°' r - 2 > and k a (5 t ) = 10 • S"" 1 

together with 

«A(P r ) = P 1 " 1 k a(p) 
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for all remaining plateau- free primes. It has been conjectured that this covers all primes [44] . 
No exception is known to date; the conjecture was tested for all p < 10 8 in [5j. Note that 
each individual prime can be analysed on the basis of Proposition CD 
The periods mod p are k a (2) = 3, k a (5) = 10, together with 



k a(p) = 




for odd primes p ^ 5, where (|) denotes the Legendre symbol and m p € N is a characteristic 
integer that covers the possible order reduction. It is 1 in 'most' cases (in the sense of a 
density definition), but there are infinitely many cases with m p > 1; this integer is tabulated 
to some extent in [441 [25] . 

Let us write down the generating polynomials for the distribution of cycles on the lattices 
L n . Once again, this is only necessary for n a prime power. We use a formulation with a 
factorisation that shows the structure of orbits on L p r \ L p r-i. In the notation of [12] . one 
finds Z x {t) = (1-i) and 

Z 2 ,(t) = (l-t)(l-t 3 )n(l-t 3 - 2 T f 

1=0 

with r > 1 for the prime p = 2, as well as 

^^t) = (i-t)n((i-^)(i-t io - 5f )) 2 ^ 

1=0 

with r > 1 for p = 5. As usual, we adopt the convention to treat an empty product as 1. The 
remaining polynomials read 

y 1 2 

Z p r(t) = (1-t) Y[(l-t K A^P e )^ pi , 

e=o 

as long as the plateau phenomenon is absent (see above). 

A. 2. Fibonacci cat map. Closely related is the matrix Mp = (JJ) e GL(2,Z), which 
is the unique square root of the Arnold cat map Ma in GL(2,Z). It appears in numerous 
applications; see |38] EJ [TJ [16] and references therein for some of them. Here, the iteration 
numbers are the Fibonacci numbers themselves, and the periods are the so-called Pisano 
periods; compare [42, A001175] and references given there, or [44] . 

The matrix Mp is not reversible in GL(2,Z) (while its square M\ is, see above), and has 
the same symmetry group as Ma- In fact, ±Mp are the only roots of Ma in GL(2,Z). This 
situation implies that the orbit structure for Mp must be such that the iteration of its square 
gives back the counts we saw in the previous example. 

For prime powers p r , with r G N, one finds K F (5 r ) = 20 • 5 r_1 together with 
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for all remaining primes, with the same proviso as for the Arnold cat map. The periods n F (p) 
are given by k f (2) = k a (2) = 3 together with 

K F (p) = 2k a (p) 

for all odd primes, which is not surprising in view of the relation between the two matrices 
M F and M A . 

The orbit distribution is more complicated in this case, as usually orbits of two possible 
lengths arise in each step. One finds 

z^t) = (i-i)^!-^) 2 ' 

and 

for the primes 2 and 5 (with r € No as before), as well as 

y> \ 2 

Zpr (t) = (i-t) J](i-a-F(pV) 2 ^(i-^F(pV)^y p '" np 
e=o 

for all remaining primes that are free of the plateau phenomenon (which possibly means all, 
see above). Here, n p € No is a characteristic integer which often takes the values 1 or 0, but 
does not seem to be bounded. 
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